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ABSTRACT 

Context. The standard dynamo model for the solar and stellar magnetic fields is based on the aCl mechanism, namely, an interplay 
between differential rotation (the CI effect) and a mean electromotive force generated by helical turbulent convection flows (the a 
effect). There are, however, a number of problems with the a effect and aCl dynamo models. Two of them are that, in the case of the 
Sun, the obtained cycle periods are too short and the magnetic activity is not sufficiently concentrated at low latitudes. 
Aims. We explore the role of turbulent induction effects that may appear in addition to the a effect. The additional effects result from 
the combined action of rotation and an inhomogeneity of the large-scale magnetic field. The best known of them is the Clx J effect. 
We also include anisotropic diffusion and a new dynamo term which is of third order in the rotation vector CI. 

Methods. We study axisymmetric mean-field dynamo models containing differential rotation, the a effect and the additional turbulent 
induction effects. The model calculations are carried out using the rotation profile of the Sun as obtained from helioseismic measure- 
ments and radial profiles of other quantities according to a standard model of the solar interior. In addition, we consider a dynamo 
model for a full sphere which is solely based on the joint induction effects of rotation and an inhomogeneity of the large-scale mag- 
netic field, without differential rotation and the a effect (a S 2 dynamo model). This kind of dynamo model may be relevant for fully 
convective stars. 

Results. With respect to the solar dynamo, the inclusion of the additional turbulent induction effects increases the period of the dy- 
namo and brings the large-scale toroidal field closer to the equator, thus improving the agreement of the models with the observations. 
For the S 2 dynamo working in a full sphere, we find dynamo modes which are steady if the effect of anisotropic diffusion is not 
included. The inclusion of anisotropic diffusion yields a magnetic field oscillating with a period of the order of the turbulent magnetic 
diffusion time. 
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1. Introduction 

Most solar and stellar dynamo models use the scenario proposed 
by Parker ( |1955| |1979| l where the magnetic field is produced by 
an interplay between differential rotation (the CI effect) and the 
collective action of turbulent cyclonic convection flows, widely 
known as the a effect (Steen beck et al.|[l966 l Krause & Radler 
1980). The scheme suggests that the a effect is responsible for 
the generation of the poloidal component of the large-scale mag- 
netic field (LSMF) of stars and other cosmic bodies. 

There is an ongoing debate on a number of problems 
connected with the a effect and aCl dynamos (see, e.g., 



Ossendrijver 2003 ; Rudiger & Hollerbach 2004, Brandenburg 
& Subramanian |2005[ ). For instance, the period of the solar 
activity cycle poses a problem. Namely, for mixing-length es- 
timates of the turbulent magnetic diffusivity in the convection 
zone and dynamo action distributed over the whole convection 
zone, the obtained cycle periods are generally much shorter than 
the observed 22 yr period of the activity cycle. For thin-layer dy- 
namos the situation becomes even worse. 

Furthermore, according to standard theory the strength of the 
a effect follows roughly a cos 6 colatitude profile. That is, the 
effect is strongest near the poles. Similarly, also the variation of 
the solar rotation rate with radius, responsible for the CI effect, 
is strongest at high latitudes ( Schou et al. 1998 ). But solar mag- 
netic activity in the form of active regions is mainly observed in 



latitudinal belts relatively close to the equator. As a possibility 
to bypass this discrepancy, meridional (poloidal) flows are un- 
der discussion, leading to so-called flux-transport dynamos (see, 
e.g., |Ossendrijver| 12003 1 Rudiger & Hollerbach| |2TJ04 Dikpati 



& Gilman, 2007jk Such flows may transport toroidal magnetic 



flux toward the equator and their speed may determine the cycle 
period. 

In this paper, we consider the possible role of a turbulent dy- 
namo mechnisms that may complement the a effect or may be 
an alternative to it. Namely, according to mean-field dynamo the- 
ory, there are other turbulent sources of the LSMF besides the a 
effect in rotating electrically conducting fluids. The influence of 
the turbulence on the LSMF is expressed by the mean turbulent 
electromotive force (MEMF) £ = (u x b), where u and b are 
the fluctuating parts of the velocity and magnetic field (angular 
brackets denote averages). Here we investigate the role of the 
Six J effect (Radler 1969 1 in axisymmetric mean-field dynamo 
models (£1 is the angular velocity of the stellar rotation and J 
the large-scale (or mean) electric-current density). This effect, 
which may be interpreted as resulting from an anisotropic tur- 
bulent electrical conductivity, has been little investigated in the 
context of solar and stellar dynamos. In the commonly used rep- 
resentation of the MEMF on the basis of symmetry arguments 
(see|Radler|[T980l|Krause & Radler|[T9"80"l|Radler| [20001 |Radler| 
|etal.||2003 



, the Six J effect represents a contribution to the 6 
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term; another contribution to this term is the shear-current effect 



( |Rogachevskii & Kleeorin||2003j|2004 i. 

An illustration of the SI X J effect is given in Fig. [T] The 
vector of the global rotation points in the vertical direction, and 
the horizontal field lines of a purely toroidal LSMF (solid blue 
lines) are linked by loop-like field lines (solid and dashed black 
lines) of a magnetic perturbation. The strength of the toroidal 
LSMF (represented by the field line density) decreases in the 
direction of the rotation vector. The Six J effect comes about 
via the following steps. First, there is a contribution (b ■ V)(B) 
to the small-scale Lorentz force, giving rise to a velocity pertur- 
bation u ~ (b-V) (B) parallel to the toroidal LSMF, indicated 
by dashed red arrows. Second, the Coriolis force deflects the 
velocity perturbation in the direction perpendicular to SI and 
(B), that is, in the direction perpendicular the plane of the fig- 
ure in the schematic (in the radial and latitudinal directions on 
the Sun or star). Third, the resulting deflected small-scale flow 
and the small-scale magnetic field give rise to a mean electro- 
motive force £ ~ <[(*-V) (B) xSl]xb) ~ - (b 2 ) (SI • V) (B), 
where in the last step we have assumed that b is mainly parallel 
to £2, so that a component of £ parallel to SI could be neglected. 
The induced electromotive force is then parallel to the LSMF. 

The MEMF was calculated analytically within a sim- 
plified version of the 
Subraman ianj |2005 



approximation (cf. 
by [Prpiiil (|2008>. The 



Brandenburg & 



relevant part of 

the MEMF, containing the combined effects of rotation and a 
nonuniformity of the LSMF, reads 



+ f { f ) s inm e n e l (2sBi m - (e + 1) B m ij 



(1) 



+ ejfe n B„ 



« (0)2 W, 



where e = fi/Jfi| is the unit vector in the direction of the ro- 
tation vector, Bjj = d{Bj)/dxj the gradient tensor of the mean 
magnetic field (we here use Cartesian coordinates xi,X2,x^ and 
the summation convention), m (0) the small-scale or turbulent con- 
vective background velocity as present in the absence of rota- 
tion and a mean magnetic field, t c the correlation time of « (0) , 

s = ^(b (l>>2 ^/ (u c V/Uop) the square root of the ratio between 

the energies of a fluctuating magnetic background field # (0) , as- 
sumed to be generated by a small-scale dynamo, and the back- 
ground velocity field « (ff) (u c - ^j(u^ 0)2 ^ is the rms value of the 

latter one and p the mass density), and and fr* (i = 1,2,...) 
denote functions of s and the Coriolis number LT that are given 
in the Appendix. Q*, measuring the influence of rotation on the 
turbulence, is defined by ST = 2Qot c , with SIq denoting the 
solid-body rotation rate; in our numerical calculations for the 
Sun we have chosen the value of the equatorial angular veloc- 
ity at the solar surface for Qq. Differential rotation, which gives 
rise to the shear-current effect, is not included at this place. We 
assume energy equipartition between the two background fields, 
i.e., e = 1. 

By its nature, the assumed small-scale dynamo is fully inde- 
pendent of the mean field. An underlying assumption here is that 
developed turbulence in an electrically conducting, non-rotating 
fluid will be magnetohydrodynamic, rather than purely kinetic, 
if the small-scale motions are sufficiently complex. With rota- 
tion added, a mean-field dynamo can be superposed to the small- 
scale dynamo. 

Hereafter, for simplicity, by the SI x J effect we mean all 
induction effects in Eq. ([1} which are of odd order in e, that is, 



the terms with coefficients ff\ ffl and respectively. Of the 
remaining two terms, with coefficients ffi and f^"\ respectively, 

that with coefficient f!f* corresponds to the B term in commonly 
used representations of £ and describes an isotropic turbulent 
diffusion. The action of the term with coefficient f[ a) will be re- 
ferred to as anisotropic diffusion. The second part of this term, 
proportional to /f a) (e+ 1), is known to describe an extra diffusion 



along e ( Kitchat movetaLl [1994} [Kitchatinov] [2002] |2004} . 

The Six J effect as it is usually understood is contained in 
the two source terms of first order in e, with coefficients /j and 

, respectively. Namely, in vector form the sum of these two 
terms can also be written as 

fl d) V(e-(B)) + sf^(e-V)(B) 

= ffV(e ■ (B)) + sjf[V(e ■ <B» - e x (V x <B»] (2) 
= -effe x m J + (jf + ejf) V(e ■ <B» . 

By their form, the two terms on the right-hand side of the last 
of Eqs. (|2]) correspond to the <5i (here, Six J) and 62 effects, re- 
spectively (cf. Krause & Radler 1980 1. In the axisymmetric case 
considered here, where gradients in the azimuthal direction van- 
ish, the 62 effect does not contribute to the important azimuthal 
component of the MEMF. Furthermore, it becomes a gradient, 
and its mean-field induction effect thus vanishes completely, if 
the Coriolis number ST (i.e., the coefficient of V(e ■ (B))) does 
not vary spatially (as in a model considered in Sect. 3.2 below). 

The term of third order in e in Eq. |lj, with coefficient f^ d \ 
is an additional source term that does not seem to have been 
included in dynamo studies before. It involves only the symmet- 
ric part of the gradient tensor of the mean magnetic field (since 
e„e m B mn = e n e m [b,,,,, + B nm j /2) and the coefficient of this sym- 
metric part, e/ 3 (£/) e,e„e m , is symmetric in the indices m, n. By its 
formal structure it thus belongs to the k term in respresentations 
of the MEMF as used in more recent papers of Radler and his 
collaborators (see, e.g., 



Radler 2000 Radler et al. 2003 Radler 



|& Stepanov|[2006l l. 

First results on the Six J effect in mean-field dynamo mod- 
els were, for instance, given in |Radler| ( |1969| l, |Stix| ( |1976[ ) and 
[Krause & Radler ( 1980). At that time the dependence of the solar 
rotation rate on radius and latitude was not known well enough. 
In the early papers also the possibility of combining the aSl and 
6S1 mechanisms was discussed. Recently, the idea of a combi- 
nation of the a and Six J effects in dynamo models was anew 
suggested by Kitchatinov] ( |2004] l. Here we explore axisymmet- 
ric kinematic a z 6Sl dynamo models (a effect plus Six J effect 
plus differential rotation) for a convective spherical shell and for 
a full sphere, following the suggestion of Kitchatinov ( 2004J and 
using the calculations of Pipin (2008 1. 

In the context of the model the following points concerning 
the Six J effect are important: 1) The effect generates a MEMF 
along the LSMF, similar to the a effect. 2) The strength of the 
Six J effect depends both on the intensity of the fluctuating mag- 
netic fields and on their linkage with the LSMF. The latter con- 
tributes to the amount of magnetic and current helicities in the 
volume considered. Thus, the strength of the Six J effect, as well 
as that of the a effect, is connected with the evolution of these 
quantities. In a companion study (Pipin, 2007), the nonlinear sat- 
uration of the a and Six J effects was investigated by integrating 
the mean-field equations coupled to an evolution equation for the 
small-scale current helicity forward in time. Similar to the non- 
linear back-reaction of the mean magnetic field on the a effect, 
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known as a quenching, the fix J effect may be suppressed by 
strong mean fields. Here we concentrate on the linear-stability 
problem. Rather than numerically simulating time evolutions, 
the eigenvalues of the Jacobian matrix of the equations for the 
mean magnetic field are calculated directly. For a comparison 
of model calculations for the a effect and helicities in the solar 
convection zone, carried out using similar approximations as for 
the calculation of £ in the repesent study, with observations at 



atmospheric levels, we refer to Kuzan yan et al.| (f2006) 

We construct an a 2 6Ll dynamo model for the Sun with dis- 
tributed dynamo action in the bulk of the convection zone. Such 
distributed dynamos have to be distinguished from boundary- 
layer dynamos that operate in the overshoot layer at the bot- 
tom of the convection zone, a thin transition region between 
the convection zone and the convectively stable radiative core, 
which is believed to coincide with the tachocline, where the 
differential rotation changes into rigid rotation in the radia- 
tive core. The main argument in favour of boundary-layer dy- 
namos is that in the convection zone proper, due to the action 
of magnetic bouyancy, magnetic flux might not be stored long 
enough to allow the generation of a sufficiently strong toroidal 
field by differential rotation. This argument is based on the pic- 
ture that the magnetic flux is concentrated in thin flux tubes 
(see, e.g. [Spiegel & Weiss) [1980} |Schiissler| [1980) |Galloway 
& Weissj 1981) |Schussler & Ferriz-Masj |2003|). Large-scale 



magnetic fields are only little affected by magnetic bouyancy 
(Kichatinov & Pipin 1993] >. A critical discussion of arguments 
for and against deep-seated and distributed dynamos, respec- 
tively, is found in Brandenburg (2005 ). 

In addition, we develop a kinematic model of an axisym- 
metric 6 1 dynamo in a full sphere. This dynamo is solely based 
on the turbulent electromotive force £> id) given by Eq. (JTJ. It 
might be working in fully convective stars. These objects pre- 
sumably do not possess layers with a strong velocity shear like 
the tachoclines at the bottom of the convection zones of solar- 
type stars, where the toroidal part of the large-scale magnetic 
field is believed to be generated. Nevertheless, fully convective 
stars are scarcely less magnetically active then stars with a radia- 
tive core and should, thus, harbor dynamo action. Discussions 



of the resulting dynamo problem may be found in the recent 
studies of |Dobler et al.| ( |200"6l >, |Chabrier & Kiiker| ( [2006> and 
Browning (2008|. Mean-field dynamo models for fully c onvec 
tive stars have been mainly based on the a 1 mechanism ( Kiiker 



& Riidiger| [1999) |Elstner & Riidiger| [2007) |Chabrier& 



Ciiker 



2006 ). Also here, we suggest the fix J effect (in its extendec 
form) as a complement or an alternative to the a effect. 

The remainder of the paper is organized as follows: In Sect. [2] 
we describe our models of the solar a 2 Sfl dynamo and the 5 2 dy- 
namo, as well as the used numerical procedure. Then, in Sect. [3] 
we present the obtained results. Finally, in Sect. [4] we draw con- 
clusions and discuss our results. 



2. Models and numerical procedure 

2.1. The models 

We restrict ourselves to axisymmetric models. Mean fields may 
well be defined by azimuthal averages. Also, the large-scale so- 
lar magnetic field is largely axisymmetric; the deviations from 
axisymmetry, as they appear, for instance, in the form of active 
longitudes, are in general small. Nevertheless, a more compre- 
hensive study would have to check whether our model is stable 
against non-axisymmetric perturbations. 

The axisymmetric LSMF is represented in the form 



(B) = curl 



r sinfl 



+ Be* 



(3) 



as the sum of a poloidal and a toroidal component; A(r, 6, t) and 
B(r, 9, t) are scalar functions of radius r, colatitude 9 and time 
f, and #0 is the unit vector in the direction of the azimuthal co- 
ordinate (p. The mean-field induction equation, containing the 
effects of differential rotation, expressed by the rotation rate 
Q(r, 9) = \fl(r, 6)\, and of the MEMF, £, then takes the form 



dA 

Ik 
BB 



r sintffirf, 



_ 1 <9(Q,A) 1 
~dl ~ ~r d(r,G) + r 



dr&a 



~d~9 



(4) 
(5) 
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2.1 .1 . The solar crSQ, dynamo 



In the solar a SLl dynamo model, to keep the numerical effort 



(i.e., the spectral resolution, see Sect. 2.2 below) at a minimum, 
only the azimuthal 11 x J effect is taken into account. This may 
be justified by the fact that the toroidal part of the solar LSMF 
is much stronger than the poloidal one. However, the other parts 
of the MEMF (a effect, isotropic and anisotropic turbulent diffu- 
sion, turbulent pumping) are included in all components. Using 
the results of Pipin] ( |2008| l, the components of the MEMF in 
spherical coordinates become 



/ 2 (rf) + (l+e)/ 1 w sm z 0dsm0B 



&r = flT { ~ 



r sinfl 



80 



■ Gsin2#/| ' B 



2r 



8r 



+ C a 



G(/f cos 2 + ^ + 2/^) 
+ cos 2 + /<f + 2/f) 

+ (/rcos 2 e + / 8 (fl) -/ 9 (a) ) 



co&0 8A 
r 2 sin6» ~80 (6) 
cos# 8A 



r 2 sin6» 80 
G 8 A 

r dr 
UdA 
r dr 



/ 2 W + (1 +e)f\ u> cos 1 6 8rB 



(a). 



T] T - 



dr 



(l+s)f{" ) cos0dsm6B 
r 80 



+ G(cos 2 9 - sin 2 0)j* a) ] B 

G cos 9 8A 



r 



(7) 



sin dr 



+ (/fcos 2 . +/ r + /r)^ 

^ 4 711 I sinfl dr 

+ (/rcos 2 e +/ r +/ r)^| 



&<b — 



m 

rsinfl 



d 2 A 

'dr 2 

sin0 d 1 dA 
~r 2 ~d0sin0~d0 



+ (,-Dstn 2 .)]^,^ +j f\2s 
+ (1 -e)sin 2 6»)] 



(l-s)f[ a) (3sin20dA d 2 A 
H sin 20- 



2r d0 



drd0 



+ sin 2 6»^-| + sin 20f\ a) [(s - 1) U + sG] ^ 
dr j 1 d9 

- [/i (o) (G - (Ge + (e - 1) f/) sin 2 (?) 
+ G/r + ( £ -l)f//' fl) ]f 



(8) 



+ C„Br sin 20G/^ 



+ C 6 ef\ 



■(d) 



rsin26» 



55 



— sin v— 

dr d0 



Here G = (d/dr)logp and t/ = (d/dr) log (m 2 ) are the scale 
factors of density (p) and turbulence intensity, respectively, and 
r\T = C n t]t, with t]t = u c ( c l3, where t c denotes the correlation 
length of the background turbulence. C n < 1, C a < 1, C$ < 1 
are parameters to control the relative strengths of different tur- 
bulence effects. C a and Cg weight the a effect and the £1 x J 
effect, respectively, and the Prandtl-like number C n regulates the 
turbulence level, with which all contributions are equally scaled. 
Below, we consider the case C n = 1/5. 

A remark concerning the introduction of weighting factors 
for the a and Six J effects seems to be in order. First, we wish to 
study the dynamo onset. For the conditions of the Sun, the two 
effects have therefore to be reduced in their strength. Second, 
the t approximation, which is used to calculate £>, is based on 
heuristic closure assumptions for the turbulence; for a critical 
analysis of the approximation we refer to Radler & Rheinhardt 
(2007). Therefore, we have left a freedom to adjust at least the 
relative strengths of different turbulence effects. 

The integration domain is radially bounded by r = Q.12R Q 
and r = 0.96R Q , where the boundary conditions are A = 0, 
drB 

—- — at the bottom boundary (a usual approximation to 
dr 



perfect-conductor conditions, see, e.g., Kohler 1973 Jouve 
|et al.[ |2008| l, and vacuum conditions (that is, B = and con- 
tinuous match of the poloidal field component to an exterior po- 
tential field) at the top boundary. 

In our numerical calculations we have used a dimensionless 
form of the equations, substituting r = xR Q and t — > rjot/R^, 
where 770 is the maximum value of 777- in the convection zone. 
The estimated value for 770 is ~ 10 9 m 2 /s. The turbulent magnetic 
diffusion time on the basis of this value and the solar radius, 
our time unit, is about 15 yr. Thus, in order to match the solar 
conditions, the models should give cycle periods of the order of 
the turbulent diffusion time. 

For the construction of the model the current knowledge of 
the rotation rate in the convection zone ( |Schou et al.| |1998| l was 
taken into account. In the numerical calculations, the rotation 



profile was approximated by (cf., e.g., Godier & Rozelot 2000 1 



Q(x, 0) = Q fix, 0) 



(9) 
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with 



/ (x, 9) = ^ [435 + 50 (jc - jc ) + 220 (jc) (l - 5 cos 2 9) 



where 



- 3.5(l - 14 cos 2 + 21 cos 4 6)], 



4>{x) = 0.5 [1 + erf(50(x - x ))] 



(10) 



(11) 



and xq = 0.71 is the position of tachocline, situated below the 
bottom boundary of the integration domain. 

The radial profiles of characteristic quantities of the turbu- 
lence, such as the rms value u c and the correlation length and 
time l c and t c of the convective background velocity field, as 
well as the density stratification were calculated on the basis of 
a standard model of the solar interior (Stix, 2002). 

Some basic model quantities, namely, the differential rota- 
tion, the turbulent pumping velocity for the toroidal component 
of the LSMF and the radial profiles of the Coriolis number Q* 
and of sec 9 (a^ measures the strength of the azimuthal 
a effect) are shown in Fig. [2] (see also Seehafer et al. 2003 
Kuzanyan et al. 2006| l. 



2.1.2. The 6 1 dynamo 

In our 8 2 dynamo model for a full sphere, representing a star, dif- 
ferential rotation and the a effect are neglected. For simplicity, 
we assume a constant density stratification R+G = -20 through- 
out the star (/?* is the radius of the star), corresponding to the 
value of G at a distance of R*/2 from the center of an M5 
dwarf with mass M - Q.2M Q , luminosity L = 0.0058L©, radius 
R+ = Q.23R Q and surface temperature r sur f — 3458 K (we used 



the stellar evolution code TWIN (Eggleton 1971 Eggleton & 



Kiseleva-Eggleton 2002) in the frontend- version WTTS (Izzard 
& Glebbeek 2006)). We also neglect the effects of the spatial 
inhomogeneity of the turbulence and set f\j = 770, with 770 denot- 
ing the value of the homogeneous turbulent magnetic diffusivity. 
Furthermore, we set Cg = 1 . At the surface vacuum boundary 
conditions are imposed. 

Different from the treatment of the solar a 2 6Ll dynamo 
model in Sect. |2.1.1| now also the contributions of the 11 x J 
effect to the r and 9 components of £ are taken into account. The 
equations for the poloidal potential A and for Eqs. Q and 
([8]), remain unchanged (except for the omission of the term pro- 
portional to C a and U — in Eq. ([8])), while in the equation for 
the toroidal potential B, Eq. Q, the first term on the right-hand 
side, describing the effect of the differential rotation, vanishes. 
We write the equation for B in the form 



8B 
~dt 



1 



(dr(& g + &f) d{& r + &f ) y 



8r 



39 



(12) 



where S r and Sg are given by Eqs. (|6]i and (|7]i in Sect. 2.1.1 (with 



fjr — i]o, U — 0, C a — 0) and £>y and S g °' denote the additional 
contributions due to the fix J effect. For these we have 



( cos 9 (c 3 jf cos 29 + ff) I d 2 A ! M 

&r -Wei T ^— Q 

r z sin 9 \ 



{ drd9 r 89 



(13) 



(c 3 jf cos 2 + /f) 



1 d 2 A dA\\ 
~r~d& + dr)\ 



MS) 



rj s< 



(</) 



cos 



20) I 3 2 A 
\drd9 



ldA\\ 
r 89 [ 



sin 29 



2r 2 



1 8 2 A 8A 



r 89 2 



(<0. 



Se)d 2 ; 



(14) 



r sin 9 



dr 2 



C3 is a weighting factor for the component of the fix/ effect with 
coefficient (i.e., the component of third order in e). There is 
no dynamo with C3 = 1. The meaning of this coefficient for the 



model will be explained further in Sect. 3.2 below. Dynamo ac 



tion sets in when both C3 and the Coriolis number exceed thresh- 
old values. In our numerical calculations for this model, length 
and time are normalized to R+ and R\/t]q, respectively. 

2.2. The numerical procedure 

The eigenvalue problem is treated by means of a Galerkin 
method. We seek the so lutions to Eqs. @§5§ for the a 2 5Q. dy- 
namo and to Eqs. ( |4|l2[ ) for the S 2 dynamo in the form 

A (x, 9, t) = e At J] J] a nm sin 9S (x) P l m (cos 9) , (15) 

n m 

B (x, 9, t) = e*£X b nm S { n B) (x) Pj„ (cos 9) , (16) 



where S and S are linear combinations of Legendre poly- 
nomials and P) n is the associated Legendre function of degree m 
and order 1 . By these expansions the regularity of the solutions 
at the poles 9-0 and 9 — n, where both A and B are set to zero, 
is ensured. To satisfy the conditions at the radial boundaries, we 
use the "ba sis recombination" of the Legendre polynomials (see 
Boyd 2001]!. In the case of the solar a 2 6Ll dynamo this reads 

(17) 



S% M = p n to + ai^K+i to + a 2 P„ +2 (x) , 
S ( n B) (x) = P„ (x) + b x P n + x (jc) + b 2 P n+2 (x) , 
where P„{x) is the Legendre polynomial of degree n and 



h 



2n + 3 



(n + 2) 2 + 2m 

An+A 
2n 2 + 6« + 3 



a 2 - - 



b 2 



(n + l) 2 + 2m 
(n + 2) 2 + 2m 
2n 2 + 2n - 1 
2n 2 + 6n + 3 ' 



(18) 



(19) 



(20) 



For the 5 2 dynamo we exploit the symmetry of the problem, as 
(x, 9) and (-x, 9 + n) represent the same point. The boundary 
conditions on x are then satisfied with 

s$£ (*) = *(*,(*)- 2 + 2 '« + "(" + 1 > 



\ 2 + 2m + (n — 2)(n - 1) / 

S^{x)=x(P n (x)-P n - 2 (x)), (22) 

and the summations in Eqs. ( [T5|jT6| run over even m + n only 
and start from n-2. 

Integrations over radius and latitude, necessary to calculate 
the expansion coefficients a nm and b nm , were done by means of 
the Gauss-Legendre procedure. The used computer code was de- 
veloped employing the free computer algebra system Maxima, 
and the eigenvalue problem was solved with the help of Lapack 
routines, which are accessible within Maxima. For the solar dy- 
namo problem we used the first 7 modes in the radial expansion 
and the first 20 modes in the latitudinal expansion, while in the 
case of the 6 2 dynamo a decomposition with 16x16 modes was 
applied. The results were confirmed by a number of runs with 
higher resolutions. 
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RADIAL DISTANCE 

Fig. 2. Basic model quantities. Left panel: Contours of the rotation rate in the solar convection zone (left) and geometry of the 
pumping velocity for the toroidal part of the LSMF (right). Right panel: Radial profiles of the Coriolis number Ll* (solid line) and 
Q'00 sec 8 (dashed line), changes sign near the bottom of the convection zone. 



3. Results 

3. 1 . The solar a 2 5Q. dynamo 

3.1 .1 . Effects of anisotropic diffusion and the small-scale 
dynamo in the a 2 Q. dynamo 

From Eq. ([T]) or Eqs. |6]-[8]l it is seen that for the case of energy 
equipartition, s = 1, the small-scale dynamo makes significant 



contributions to different parts of the MEMF. |Kitchatinov (2002 



2004) has found that anisotropic diffusion may be in a large 
part responsible for the equatorial drift of the toroidal LSMF. 
This drift is actually an extra diffusion (of both the toroidal and 
poloidal parts of the mean field) along the rotation axis e. In 
Eq. ([T]l it is described by the second part of the anisotropic - 
diffusion term, which is proportional to f\ a \s + 1). Thus, for 
s — 1 the drift is enhanced by a factor of two compared to the 
case of s — 0. The first part of the anisotropic-diffusion term in 
Eq. jlj can in vector form be written as 2f^se x V(e • (B)). Due 
to the assumed axisymmetry, it is purely toroidal (azimuthal). 
Therefore, it does not influence the toroidal part of the mean 
field. It actually describes a diffusion of the poloidal part of the 
mean field in the direction perpendicular to the rotation vector e. 

In Fig.[3](top and middle), showing butterfly diagrams in the 
form of isocontours of the toroidal LSMF (integrated over depth 
in the convection zone) and overlaid greyscale plots of the radial 
LSMF at the top boundary, the effects of anisotropic diffusion 
and the small-scale dynamo are demonstrated for a pure a 2 Ll dy- 
namo model on the basis of the first unstable dipolar eigenmode. 
By dipolar/quadrupolar modes we mean modes that are antisym- 
metric/symmetric with respect to the equatorial plane. The used 
dipolar mode is not the primary dynamo eigenmode, cf. Fig. [4] 
(left) in Sect. 3.1.2 below. Qualitatively, the results presented 



here resemble those of Kitchatinov (2002), who did not include 
the effect of the small-scale dynamo and used a slightly different 
rotation law. 

As is seen in Fig.|2](left), the direction of the turbulent pump- 
ing velocity is such that it transports magnetic flux toward the 
equator. So also the pumping effect offers a possibility to ex- 



plain the observed latitudinal distribution and drift of the solar 
activity phenomena (Guerrero & de Gouveia Dal Pino 2008 1. In 
our model the pumping effect is obviously to weak to produce a 
noticeable equatorward drift of the toroidal LSMF. 

Finally, as is seen in Fig. [3] (bottom), the dynamo wave 
moves mainly radially outward. This is in accordance with the 
tendency of the ail dynamo waves to move along isorotation 
surfaces ( Yoshimura] |1975| l. In the bulk of the convection zone 
these surfaces are largely parallel to radius (cf. Fig. [2] left panel). 
A better agreement with the observations could thus be expected 
if the toroidal LSMF were confined to the tachocline, where 
the rotation rate varies mainly with radius. This would then 
amount to a boundary-layer dynamo, which is basically different 
from the distributed convection-zone dynamo considered here. 
A more complete model of the solar dynamo would include the 
dynamo effects of both the tachocline and the convection zone 
proper. 



3.1.2. The a 2 6Q. dynamo 

Figure |4] (left) shows a linear-stability diagram for the a 2 6Ll dy- 
namo, i.e., the stability boundary for the most unstable eigen- 
mode in the Cg-C a plane. The a 2 Ll dynamo is included as a lim- 
iting case (Cg = 0). In the figure, regions where the eigenmode 
with the largest growth rate is a dipolar mode are indicated by 
shading. In the limiting case of Cg = the dipolar modes are not 
the most unstable modes. Though a dipolar mode becomes un- 
stable only slightly above the total-stability boundary, there is a 
preference for quadrupolar modes in the pure <t 2 Q dynamo, con- 
trary to the antisymmetric parity characterizing the large-scale 
solar magnetic field. The parity issue is met in other dynamo 
models, like, e.g., the flux-transport models, as well (see, e.g., 
2005 [ l. For the pure d>Q dynamo (11 x J effect 



Dikpati et al 



plus differential rotation), on the other hand, there is a prefer- 
ence of dipolar modes (C a — in Fig. |4|. These are, however, 
non-oscillatory. An example of such a mode is shown in Fig. [4] 
(right). The fact that in the pure model only non-oscillatory 
modes are found seems to indicate that for models of the solar 
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Fig. 3. Effects of anisotropic diffusion and the small-scale dynamo for an a 2 Q. dynamo model on the basis of the first unstable 
dipolar eigenmode (C a = 0.087). Top and middle: Butterfly diagrams in the form of isocontours of the toroidal LSMF (integrated 
over depth in the convection zone) and overlaid greyscale plots of the radial LSMF at the top boundary without (top) and with 
(middle) the effects of anisotropic diffusion and the small-scale dynamo. Time is measured in units of R%Itjq. Solid/dashed lines and 
bright/dark areas indicate positive/negative field values. Bottom: Snapshots of the strength of the toroidal LSMF (greyscale plot) and 
field lines of the poloidal LSMF over half a cycle for the model containing the effects of anisotropic diffusion and the small-scale 
dynamo. Solid/dashed lines indicate clockwise/counter-clockwise field direction along the poloidal field lines. 



dynamo the a effect is also needed, at least for models without 
meridional circulation. 

For increasing C„, oscillatory dipolar modes are excited if 
C a is no longer very small compared to Cg. The properties of 
these modes are determined by the ratio C a /Cg. Fig. |5]shows a 
simulated butterfly diagram (isocontours of the toroidal LSMF 
with an overlaid greyscale plot of the radial LSMF at the top 
boundary) for Cg = 0.07, C a = 0.01. For these parameter values 
we obtain a slow dynamo wave propagating from the equator to 
the pole for both the toroidal and poloidal fields, in contrast to 



the solar observations. The period of the dynamo is considerably 
increased compared to the pure a 2 Ll dynamo shown in Fig. [5] 

In the intermediate case, when the fix J effect exceeds the a 
effect while the latter one is not very small compared to the first 
one, we find a good match with the solar case. A correspond- 
ing butterfly diagam and snapshots of the toroidal and poloidal 
fields over half a cycle for Cg = 0.1, C a = 0.05 are shown in 
Fig. [6] While the wings of the obtained dynamo waves are too 
wide, other qualitative properties, such as the directions of the 
drifts of the toroidal and poloidal parts of the magnetic field and 
their phase relation (e.g., polar reversal of B r shortly after max- 
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Fig. 4. Left: Stability boundary for the a 2 6Ll dynamo in the Cg-C a plane. The stable region (no dynamo) lies below and to the 
left of the solid line. Shading indicates dominance of dipolar modes. Right: First unstable eigenmode for the steady 6Ll dynamo. A 
greyscale plot of the toroidal LSMF and field lines of the poloidal LSMF are shown. Solid/dashed lines indicate clockwise/counter- 
clockwise field direction along the poloidal field lines. 





Fig. 6. Top: Butterfly diagram in the form of isocontours of the toroidal LSMF (integrated over depth in the convection zone) 
and overlaid greyscale plot of the radial LSMF at the top boundary for an a 2 6Ll dynamo with Cg = 0.1, C a - 0.05. Time is 
measured in units of R%/t]q. Solid/dashed lines and bright/dark areas indicate positive/negative field values. Bottom: Snapshots of 
the strength of the toroidal LSMF (greyscale plot) and field lines of the poloidal LSMF over half a cycle. Solid/dashed lines indicate 
clockwise/counter-clockwise field direction along the poloidal field lines. 



imum of at low latitudes) are in good agreement with the 
observations. The period of the dynamo is shorter than but com- 
parable with the turbulent magnetic diffusion time. In the range 
1 < C$/C a < 3 the match with the solar observations is best 
(it can be seen in Fig. [4] that in this range the dipolar modes are 
dominating at the stability boundary). 

The obtained ratio Bj/Bp of the toroidal {Bj) and poloidal 
(Bp) field strengths is about 50, the ratio of the toroidal and 
poloidal field energies thus of the order of 10 3 . This comes close 
to the ratio Bj/Bp ~ 100 estimated from measurements of the 
fields in active regions [yielding the estimate Bj ~ 200 G for a 



distributed toroidal field in the convection zone, cf. |Stix| ([2002 
Sect. 8.4. 1)] and the field strength at the solar poles (yielding the 
estimate Bp ~ 1 G for the poloidal field). 

3.2. The 6 2 dynamo 

The usual fix J effect, given by the term —efjPe x hqJ in the 
expression for the turbulent emf (cf. Eq. Q), does not contribute 
to B J. Therefore, it cannot bring energy into the mean magnetic 
field and is not capable of dynamo action when working alone. 
It can yield a dynamo, however, when acting together with dif- 
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Fig. 5. Butterfly diagram in the form of isocontours of the 
toroidal LSMF (integrated over depth in the convection zone) 
and overlaid greyscale plot of the radial LSMF at the top bound- 
ary for an a 2 6Q. dynamo with Cg = 0.07, C a - 0.01. Time is 
measured in units of R%It]q. Solid/dashed lines and bright/dark 
areas indicate positive/negative field values. 



ferential rotation. The role of the effect is then to transfer energy 
from the toroidal to the poloidal field. In our model of the 5 2 dy- 
namo differential rotation is neglected (as well as the a effect). 
Energy has thus to be provided by the additional source term in 
our model, the term of third order in e. Before studying our ax- 
isymmetric model of the 5 2 dynamo in spherical geometry, we 
consider a stronger simplified model in plane geometry, in order 
to check whether such a dynamo is in principle possible. For this 
purpose, Eq. ([T| is simplified to 

£ = - m V x (B) + ci (e ■ V) (B) + c 3 e (e ■ V) (e ■ (B)) , (23) 

where c\ and C3 are constants. This MEMF contains an isotropic 
diffusion term and two terms corresponding to the fix J effect, 
one of first order and the other of third order in e (their formal 
structures agree with those of the terms with coefficients / 4 (rf) and 

jf®, respectively, in Eq. jlji). The rotation axis coincides with the 
Z axis of a system of Cartesian coordinates x, y, z, and the mean 
magnetic field is written as 



(B) = V x [a (x, z) e v ) + b (x, z) e y 



(24) 



e y is the unit vector in the y direction, in which the system is in- 
variant (corresponding to the <p direction in the spherical model). 
From the mean-field induction equation 



d{B) 
dt 



Vxfi 



we then get 



d 




dt 


-m 


a 




m 


-m 



dx 2 dz 2 



dx 2 dz 2 



db 



a = c\ 



dz ' 

b = -Cl—r- (Cl + C 3 ) 

OZ i 



d 3 a 
dx 2 dz 



(25) 

(26) 
(27) 



Solutions to Eqs. ( 26 ■ 27 1 may be sought in the form 



a — a exp {yt + i (kx + qz)} , b — b exp {yt + i (kx + qz)} , 

(28) 

with y being determined by the condition that the determinant of 
the resulting system of algebraic equations vanishes, that is, 

[y + 770 (k 2 + q 2 )f + q 2 [c\ (k 2 + q 2 ) + Cl c 3 k 2 ] = , (29) 



Fig. 7. Stability diagrams for the 6 2 dynamo in spherical geom- 
etry without the effects of anisotropic diffusion and turbulent 
pumping for C3 = -2.5 (left) and C3 = 2.5 (right). The four 
largest real parts of the eigenvalues A are shown as functions of 
the Coriolis number. Only dipolar modes are considered. 



with roots 



Tl,2 = ~m (k 2 + q 2 ) + ^-q 2 [c 2 (k 2 + q 2 ) + C1C3& 2 ] . 



(30) 



It is seen that y is real for growing modes, so that only non- 
oscillatory dynamo modes can be excited. Furthermore, neces- 
sary, but not sufficient conditions for the existence of a dynamo 
mode are (i) c x + 0, (ii) c 3 + 0, (iii) k + 0, (iv) q + and (v) 



c 3 k 1 + 

cl > ~. k 2 



(31) 



That is, ci and C3 must have opposite signs and the absolute value 
of C3 must be larger than that of C\. Now in the limit of fast 
rotation, Q* » 1 , the signs of / 3 (flf) and / 4 (rf) are opposite, namely, 

for LT — > 00 one gets ji^/ji — > -2. Suppose ci = — C3/2. Then 
the (necessary and sufficient) condition for dynamo instability 
becomes 

c\q 2 (k 2 - q 2 ) 



> 1 



(32) 



It can be satisfied if \k\ > \q\. The optimal wave number ra- 
tio is \k\l\q\ — V3, giving dynamo excitation for c 2 > 32tjI. 
These findings are largely confirmed by the kinematic 5 2 dy- 
namo model in spherical geometry, for which we now discuss 
the eigenmodes of dipolar type. 

In the spherical model we have left only one of the above 
two parameters, C3, as a weighting factor for the term with coef- 

14 1. Besides 
lfwe neglect 
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ficient (term of third order in e), see Eqs 
that we now also vary the Coriolis number £ 
the effects of anisotropic turbulent diffusion (but not those of 
isotropic turbulent diffusion) and density stratification (i.e., tur- 
bulent pumping) in the spherical model, we get steady dynamo 
modes as obtained in the plane model above. Stability diagrams 
for this case are shown in Fig. [7] where the left panel shows that 
there is no dynamo effect if the coefficients of the terms describ- 
ing the fix J effect do not have signs in accordance with what 
is expected to be necessary for dynamo action from the con 



sideration of the plane model. Namely, /- 

Ad), Ad) 



Ul) 



is negative and f. 



id) 



positive for all Q*, \ff'/f^ a '\ monotonically increasing with Ll* 

and / 3 ~ -2 for Ll* > 5. For C3 = 1 we found no dynamo 

instability in the whole £2* interval. The condition for dynamo 
action was found to be C3 > 2. Furthermore, a threshold value 
of Q* must be exceeded. The larger C3, the smaller the threshold 
value of Q.*. 

Including anisotropic turbulent diffusion decreases the 
threshold values of C3 and £2* for instability, while including 



10 



V. V. Pipin and N. Seehafer: Stellar dynamos with fix J effect 






Fig. 8. Snapshots of the magnetic field in the time-periodic 6 2 dynamo (first unstable mode) over half a cycle for C3 = 2.5, LT = 9.1 
in the form of greyscale plots for the toroidal LSMF (bright/dark areas indicate positive/negative field values) and field lines of the 
poloidal LSMF (solid/dashed lines correspond to clockwise/counter-clockwise field direction). 



turbulent transport (pumping), associated with density strati- 
fication, raises them. Yet, including one of them or both to- 
gether changes the character of the dynamo bifurcation from a 
steady-state to a Hopf bifurcation. Namely, the two most un- 
stable modes then merge, yielding a slowly oscillating dynamo 
mode. "Slowly" here means that the time period of the mode is 
larger than the turbulent diffusion time of the system. For the pa- 
rameter choice C3 = 2.5, R+G = -20, the threshold value of the 
Coriolis number is Q* » 9.1, and the frequency of the obtained 
dynamo mode is co » 0.8 inverse diffusion times. Fig. [8] shows 
the evolution of the large-scale toroidal and poloidal field com- 
ponents over half a cycle. It is seen that the strongest toroidal and 
poloidal fields are concentrated close to the rotation axis in the 
polar regions (actually the field lines of the poloidal field are iso- 
contours of the potential A, so the strength of the poloidal field is 
also seen). The dynamo oscillation closely resembles a standing 
wave. 

4. Conclusions 

We have studied kinematic axisymmetric mean-field dynamo 
models in which the 11 x J effect was included as one of the 
mechanisms by which the large-scale magnetic field is gener- 
ated. The effect was used in an extended form, together with 
all other effects that result from the combined action of rotation 
and a nonuniformity of the large-scale magnetic field, including 
anisotropic turbulent diffusion. Concurrently to the a effect, the 
II x J effect may, in particular, generate the large-scale poloidal 
magnetic field of stars. Our results show that the inclusion of 
the additional effects alleviates some of the problems connected 
with the a effect and aQ. dynamos. Compared to the standard aQ. 
dynamo, the inclusion of the fix J effect increases the time pe- 
riod of the dynamo and brings, thus, the models in better agree- 
ment with the solar observations. Furthermore, the large-scale 
toroidal field comes closer to the equator, bringing the models in 
better agreement with the observations also in this respect. The 
observed phase relation between the toroidal and poloidal field 
components is correctly reproduced and, in contrast to the aLl 
model, dynamo modes with the correct (dipolar) parity become 
unstable first. 

We did not explore the possible effects of meridional circu- 
lation in the considered models. Circulation-dominated dynamo 
models are very popular presently. They apparently avoid the 
problems with aQ. dynamos mentioned above. But they also 
introduce (yet) unproven assumptions concerning the merid- 
ional flow and work only if the turbulent magnetic diffusiv- 
ity is reduced by about two orders of magnitude compared to 
the mixing-length estimates. We do not claim that the mod- 
els considered in the present paper are superior to circulation- 
dominated dynamo models. Possibly the effects studied here 



and meridional circulation should be combined. We have pre- 
liminary results (not published yet) which indicate that a 5Q. 
dynamo model with meridional circulation, but without the a 
effect, admits oscillatory eigenmodes that give a solar-type dy- 
namo. Another possible improvement of our model would be the 
inclusion of the tachocline. 

We also gave a first example of a 5 2 dynamo which is solely 
based on the joint induction effects of rotation and an inhomo- 
geneity of the large-scale magnetic field, without differential ro- 
tation, the a effect and turbulent pumping. This kind of dynamo 
may be relevant for fully convective stars. We found a dynamo 
instability both for a strongly simplified model in plane geome- 
try and for an axisymmetric model for a full sphere. The dynamo 
modes are steady if the effect of anisotropic diffusion is not in- 
cluded. In the case of the sphere, the inclusion of anisotropic dif- 
fusion yields a slowly oscillating magnetic field. Including tur- 
bulent pumping (due to density stratification, not part of a pure 
6 2 dynamo) leads to the same result. The period of the oscilla- 
tion is of the order of the turbulent diffusion time of the system. 
For M dwarfs this is estimated to be 10 yr - 100 yr. Preliminary 
results for combinations of a 2 and S 2 dynamo models (not con- 
tained in the paper and not yet published elsewhere) indicate that 
the mixture of the two effects generally produces oscillatory dy- 
namo modes. Our aim here was not to provide a dynamo model 
for fully convective stars, but merely to suggest the combined 
effects of rotation and an inhomogeneity of the large-scale mag- 
netic field as ingredients in future dynamo models. 
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Appendix 



Here we give the definitions of the functions and ff" 1 that are 
used in the representation of the turbulent electromotive force £>. 
For details of the calculations we refer to Pipin (2008). 
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